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Transition operator method is proposed for description of the dynamics of spectroscopic transi- 
tions. Quantum-mechanical analogue of Landau-Lifshitz equation has been derived for the system 
representing itself the periodical ferroelectrically (ferromagnetically) ordered chain of N equivalent 
elements, interacting with external oscillating electromagnetic field. Landau-Lifshitz equation was 
represented in Lorentz invariant form by using Hilbert space over the ring of quaternions. It has 
been shown, that spin vector can be considered to be quaternion vector of the state of the system 
studied. From comparison with experiment for the first time from pure optical measurements the 
value of spin S — 1/2 for optically active centers - spin-Peierls solitons in carbon chains - has been 
obtained. The ratio of imagine to real components of complex charge is evaluated for given centers 
to be § « (1.1 - 1.3)10 2 . 
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I. INTRODUCTION 



Spectroscopic description of the systems, interacting 
with electromagnetic (EM) field including the descrip- 
tion of transitional effects, for instance, Rabi-oscillations, 
free induction, spin echo in magnetic resonance spec- 
troscopy and corresponding optical analogues in optical 
spectroscopy is achieved within the framework of gyro- 
scopic model. Mathematical base for gyroscopic model is 
Landau-Lifshitz (L-L) equation. L-L equation was pos- 
tulated by Landau and Lifshitz for macroscopic classical 
description of the motion of magnetization vector in fer- 
romagnets in 1935 [l]. L-L equation is: 



f = [7«Hx4 



(1) 



where S is magnetic moment, H is effective magnetic 
field, 7 H is gyromagnetic ratio. Clasical gyroscopic model 
for transition dynamics and for description of transitional 
effects was introduced formally Q by F.Bloch, who has 
added the phenomenological relaxation term to L-L equa- 
tion. Then L-L equation and Bloch equations, based 
on it, were substantiated quantum- mechanically in mag- 
netic resonance theory, but the only for describtion of 
the motion of magnetic moment vector in external mag- 
netic field 0,3- Quantum- mechanical Bloch equations 
were in fact postulated for the description of dynamics 
of magnetic resonance transitions. The development of 
gyroscopic model for the characterization of optical sys- 
tems, interacting with EM-field was done on the base of 
density operator formalism, that is quantum statistically 



The optical analogue of L-L equation is: 

f = [^x 7,4 



(2) 



where 7 B was called formally gyroelectric ratio by using 
the only its mathematical analogy with gyromagnetic ra- 
tio without however corresponding physical meaning, P 
(Bloch vector), E are some axial vectors, with compo- 
nents representing themselves various physical quantities. 
So, for instance, in Q is emphasized that the only P x , 
P v and E x , E y components of vectors P, E, correspond- 
ingly, characterize the genuine electromagnetic proper- 
ties of the system, at the same time the components P z 
and E z cannot be referred to electromagnetic character- 
istics. In other words physical meaning of Bloch vector 
and E seems in existing optical transition theory to be 
unclear. In other words, the vectors P, E entering optical 
Bloch equations are in fact some mathematical abstrac- 
tions. Really, like to magnetic moment and magnetic 
field vectors in magnetic resonance Bloch equations, both 
the vectors P, E have, according to existing description, 
necessary axial symmetry, however in artificial way. For 
instance, P has two components P x , P y of polar elec- 
tric dipole moment vector, that is, both the components 
remain without any physical meaning change. At the 
same time the third component of P has entirely dif- 
ferent physical meaning and it is considered to be the 
population difference between the energy levels, see for 
example @- Mathematically the objects, which are 
like to Bloch vector, can, naturally, exist. However, even 
mathematically, it will mean, that the spaces, which pro- 
duce the set of Bloch vectors P and the set E over the 
field of real numbers, is the only affine space and it can- 
not be metrizable. In other words, the components of P, 
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E in existing interpretation cannot represent the physical 
quantity to be a measurable single whole. Physically the 
model, in which the part of components of characteris- 
tic vector possess by different symmetry (polar in given 
case) in comparison with the final symmetry (axial) of 
vector itself seems to be also incorrect. Given conclusion 
is especially evident, if to take into account, that the 
mechanism of spectroscopic transitions' dynamics itself 
has to be independent on the nature of the centers, in- 
teracting with free oscillating electromagnetic field (EM- 
field), that is, whether they possess by magnetic or elec- 
tric dipole moments, owing to independence of magnetic 
and electric components of EM-field. Given conclusion is 
in complete agreement with experiment, since dynamics 
of the spectroscopic transitions both in otical and mag- 
netic resonance areas can be described mathematically 
by the same phenomenological Bloch equations. Let us 
remember, that the analogue of Bloch vector in magnetic 
resonance spectroscopy is total magnetic moment, that 
is pure axial vector with peer all the three components. 

Note, that most general at present description of spec- 
troscopic transitions' dynamics by means of density op- 
erator formalism has a substantial disadvantage. Really, 
it seems to be understandable, that in some experimen- 
tal cases, in particular, for local single centers, which can 
be used in nanotechnology, for instance, for elements of 
binary logic, quantum informatics, quantum computing 
with very fast rates of processes of information treat- 
ment, realized on semiconductor QDs, Josephson cubits, 
spin-qubits and so on, statistical density operator formal- 
ism can be not applicable for description of transitions' 
dynamics at all. The main reason is the following. If 
the state of quantum system is mixed state, that seems 
to be necessary for functioning of nanotechnology ele- 
ments, then the event frequency of realization of definite 
pure state in the superposition can be not corresponding 
to its probability. In other words, the well known law 
of probabilities' theory - the law of large numbers has 
to be taken into account. Moreover, by measurements 
with short times, which are comparable with the resi- 
dence time in definite pure state iof the superposition, 
the amplitude of given state cannot be determined by its 
probability for any quantum system. 

New approach to a description of spectroscopic tran- 
sitions' dynamics, which allows to overcome given dif- 
ficulties, was proposed for the first time in |8|. Given 
approach is based on transition operator method. Op- 
erator equation, describing the optical transitions' dy- 
namics (for the example of simple ID-model of quantum 
system), has been obtained by using of stated method 
[8|. It has been shown, that given equation is operator 
equivalent to Landau-Lifshitz (L-L) equation [J in its 
difference-differential form, which takes usual differential 
form in continuum limit. 

The aim of the work presented is description of tran- 
sition operator method in more details and development 
of the model, proposed in 



II. SPECTROSCOPIC TRANSITIONS' MODEL 
FOR ID-CHAIN OF EXCHANGE INTERACTING 
CENTERS 

Let us consider the system representing itself the peri- 
odical ferroelectrically (ferromagnetically) ordered chain 
of N equivalent elements, interacting with external os- 
cillating electromagnetic field (EM-field). It is assumed, 
that the interaction between elements (elementary units) 
of the chain can be described by the Hamiltonian of quan- 
tum XYZ Heisenberg model in the case of a chain of 
magnetic dipoles and by corresponding optical analogue 
of given Heisenberg model in the case of a chain of elec- 
tric dipoles. We will consider for the simplicity the case 
of isotropic exchange. Each elementary unit of the chain 
will be considered to be two-level quantum system like to 
one-electron atom. Then Hamiltonians for the chain of 
electrical dipole moments and for the chain of magnetic 
dipole moments will be mathematically equivalent. The 
chain for distinctness of electrical dipole moments can be 
described by the following Hamiltonian 

H = H° +H CF +H J , (3) 

where TiP is chain Hamiltonian in the absence of the in- 
teraction with external EM-field and between structural 
elementary units of the chain themselves. We suggest, 
that given structural elementary units are optical (mag- 
netic) centers with one electric (magnetic) dipole pro cen- 
ter. Iri is Hamiltonian, describing the interaction be- 
tween classical external EM-field and atomic chain, 7i J 
is Hamiltonian of exchange interaction between above in- 
dicated centers. 

To obtain the equation, describing the transition dy- 
namics, we will use the transition operator method. 
Given method was used in atomic spectroscopy, see, for 
instance, [f|. The essence of the method, adapded for 
our task, is the following. Let {|m;)}, I = l,n, m = a,j3 
be full set of energy states of I — th optical (magnetic) 
center of the chain in the absence of the interaction both 
between the centers themselves and with EM field. Here 
\cti) is lower state and \(3i) is upper state. They are eigen- 
states oil — th item in 7i°, I = l,n. It means, that H 
is 

JV 

W = 5^£™,K)<m„|. (4) 

v—l m 

Here m = a, (3, E mv are eigenvalues of 7i°, which cor- 
respond to the states \m v ) of v — th chain center. Let 
set up in correspondence to the states \j v ), v = l,n the 
operators 

ai m ^\ Jv )( mv \, (5) 

where j = a, /?, m = a, (3. Then, it is evident, that oper- 
ators <r+ = \[3 V ) (a v \ arc the operators, which transform 
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eigenstates \a v ) into \(3 V ) and operators a~ — \a v ) (J3 V \ 
are the operators, transforming eigenstates \f3 v ) into |a„) 
for Mv, v — Let define also the operators cr^ = 

~ (a„|. The operators b% transform for Vw, 
w = l,n the sum of states \j3 v ) + |a„) into their differ- 
ence and the reverse is also true. The above defined set of 
spectroscopic transition operators can be completed by 
unit operators , v = l,n, which, owing to complete- 
ness of state sets can be represented by the following 
relationship 



E 



m v ) (m v 



m = a,f3, and by zeroth operators crj], v 
can be defined in the following way 



<7° = al 



(0) 



1, n, which 



(7) 



where j = — , +, z, E. It is easy to show, that at fixed v 
the set of transition operators is closed relatively the al- 
gebraic operations of commutation and anticommutation 
(that is relatively taking of commutator and anticommu- 
tator) and relatively the operation of Hcrmitian conju- 
gation. For instance, the relationships for commutation 
rules are 



[d* n ,d**] = (%s mp -d*rs, 



that is, 



2a: 



2(7.1 



(8) 



(9) 



So, at fixed v the set of transition operators {<5^ m } pro- 
duces algebra. In particular, the set of {<r™} operators, 
where m is z, +, -, 0, E, produces algebra, which is iso- 
morphic to S = 1/2 Pauli matrix algebra, completed by 
unit and zeroth matrices, that is, mappings 



fv 



ap,m = z,+,-,0,25 



(10) 



realize isomorphism. Here v is a number of chain unit, 
<7p is extended set of Pauli matrices for the spin of 1/2. 
unit. 

The structure of transition operators indicates on 
mathematical advantages of transition operators' method 
in comparison with the method of density matrix. Re- 
ally, the method of density matrix gives the relationship 
for probabilities of the states, in spectroscopic transition 
operators' method we are dealing immediately with the 
states. 

Let us define the vector operator: 



(11) 



Consequently, taking into account (|10[) . from physical 
point of view a v represents itself some vector operator, 
which is proportional to vector operator of the spin of 
v — th chain unit. 

It is evident, that the set of vector operators pro- 
duce linear space over field of complex numbers, which 



can be called transition space. It is 3-dimensional (in 
the case of two- level systems), that is 3 matter oper- 
ator equations of the motion for components of a v is 
necessary for correct description of optical transitions in 
semiclassical approach. Typical inaccuracy in many of 
the quantum optics calculations of transition dynamics is 
connected with obliteration of correct dimensionality of 
transition space. Further, isomorphism of Pauli matrix 
algebra and spectroscopic transition operators' algebra 
allows to change the spin variables in known Heisenberg 
exchange Hamiltonian by components of vector a v . 

We will use further the rotating wave approximation 
[5J. Then the chain for distinctness of electrical dipole 
moments can be described by the following Hamiltonian: 



JV 



N 



■p^W+e-^ + fr-e***) 



N 1 



(12) 



n=l 



It is suggested in the model, that \a n ) and \(3 n ) are eigen- 
states, producing the full set for each of TV elements. It 
is evident, that given assumption can be realized strictly 
the only by the absence of the interaction between the 
elements. At the same time proposed model will rather 
well describe the real case, if the interaction energy of ad- 
jacent elements is much less of the energy of the splitting 
huj = £p — £ a between the energy levels, corresponding 
to the states \a n ) and \(3 n ). This case includes in fact 
all known experimental situations. Further, p"" is ma- 
trix element of dipole transitions between the states \a n ) 
and \j3 n ). Matrix elements along with energy difference 
of given states £p — £ a are suggested to be independent 
on n, Ei is amplitude of electric component of exter- 
nal EM-field on the n-th element site, h is Planck's con- 
stant, Je is optical analogue of the exchange interaction 
constant. Here, in correspondence with the suggestion, 
J„ = JZ = J l i — J 7 '- In the case of the chain of magnetic 

& E E E ° 

dipole moments electrical characteristics of the systems 
studied in Hamiltonian (|12[) are replaced by correspond- 
ing magnetic characteristics 

25 x n 

n a(3 _^ n af3 (13) 



■H[\ 



-> -g H (3 H H a = j H H , 

where 22" is amplitude of magnetic component of exter- 
nal EM-field (microwave field) on the n-th element site, 
J H is the constant of magnetic spin exchange interaction, 
matrix element of magnetic resonance dipole tran- 
sitions, Hq is external static magnetic field, /3 H is Bohr 
magneton, g H is g-tensor, which is assumed for the sim- 
plicity to be isotropic. 

The first term in Hamiltonian (fT2")) characterizes the 
total energy of all chain elements in the absence of exter- 
nal field and in the absence of interaction between chain 
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elements. The second item characterizes an interaction of 
a chain with an external oscillating EM-ficld in dipole ap- 
proximation. Matrix elements of full dipole up and down 
transitions between couples of the states (\a n ), |/3 n )) and 
(|/3 n ), |a„)), respectively, are assumed to be equal, that 
is, spontaneous emission is not taken into consideration. 
The third item is, in essence, Hamiltonian of quantum 
Heisenberg XXX-model in the case of magnetic version 
and its electrical analogue in electric version of the model 
proposed. The equation of the motion for a k is: 



ift 



dak 
dt 



[aj-,H] e + + [&+ H] e_ + [a z k ,H] e z . (14) 



We will consider the case of homogeneous excitation 
of the chain, that is E™ is independent on unit num- 
ber n {Ei = E\). Then, using the commutation rela- 
tions [a£,&~] = 5 kn a^, [6--,o-£] = 25 kn a-, [a z k ,a+] = 
ZSkn&ni we obtain 



dt 



(15a) 



2iJ 



dt 



(15b) 



da7 



dt 



-iuj <J k 



iQ„e 



(15c) 



i J 



-f (to , (H+x + + {°i > + *fc-i)» . 

where expressions in braces { , } are anticommutators. 
Here g and j E are Rabi frequency and gyroelectric ratio 
respectively. They are determined, correspondingly, by 
relationships 



a/3 



Exrff p' 



which are replaced by relationships 
Hypf 



„Q/9 



n. 



= 7„#i, 7„ = -g- 



(16) 



(17) 



in the case of the chain of magnetic dipole moments. 
The equations (|15p can be represented in compact vector 
form, at that the most simple expression is obtained by 
using of the basis e+ = (e x + ie y ), e_ = (e x — ie y ),e z . 
So, we have 



dak 
dt 



<7fc x <& h . 



l,k+l 



(18) 



where a k is given by (jlljl . but with the components in 
matrix representation of operators, corresponding to new 

basis, k = 2,7V— 1, and vector operator &k-i fe+i = 



®k-i,k+l^ + ©fc_ ljA;+1 e_ + 0|_ 1:fc+1 e z , where its com- 
ponents are 



fc+l = ^e"-* - -^Iff^ -!- rr 



2J E 



<S+ = 0, e lut - + a + 

2J 



(19) 



It should be noted, that vector product of vector oper- 
ators in (|18p is calculated in the correspondence with 
expression 



CTjfc x <5 k - 



l./v+l 



e_ x e. 



® fc-l.fc+1 



e z x e+ <r fc © fe _i, fe+ i 
e+ x e_ ^ ®Vl,fc+l 



(20) 



Given expression is similar to an expression usually used 
in vector analysis, however, anticommutators of corre- 
sponding components instead their products are taken by 
its calculation. In fact it is definition for vector product 
in the case of operator vectors. Given definition seems to 
be natural generalization of usual vector product, since 
the only in given case the result will be independent on 
a sequence of components of both the vectors in their 
products like to that one for usual vectors. Naturally, 
the expressions like to ([2"0")) can be used for calculation 
of vector product of common vectors, which is provided 
with coefficient i. 

Taking into account the physical sense of vector oper- 
ators (Tfc we conclude, that (fT8|) represent themselves re- 
quired quantum-mechanical difference-differential equa- 
tions (the time is varied continuously, the coordinates 
are varied discretely) for the description of the dynam- 
ics of the spectroscopic transitions (in the frames of the 
model proposed) . It follows from here in view of isomor- 
phism of algebras of operators a k and components of the 
spin, that the (|18[) is equivalent to L-L equation (in its 
difference-differential form) . 

Therefore, we have proved, on the one hand, the pos- 
sibility to use L-L equation for the description of the dy- 
namics of spectroscopic transitions and for the descrip- 
tion of transitional effects from unified theory positions 
for both cases, that in, in optical and magnetic resonance 
spectroscopy areas. On the other hand, we have proved 
(in the frames of the model proposed) , that the quantum 
equation for the description of the dynamics of the spec- 
troscopic transitions coincides with the equation for spin 
moment motion. Given result was used always in its im- 
plicit form by spectroscopic transition dynamics studies 
in magnetic resonance, however without its substantia- 
tion, see, for instance [3(. 

To obtain the continuous approximation of (|18p for co- 
ordinate variables, we have to suggest that the length of 
electromagnetic wave A satisfies the relationship A >> a, 
where a is lD-chain-lattice constant. It is evident, that 
in optical frequency region the relationship A >> a is 
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always takes place. Then the continuous limit will be 
realized, if to substitute all the operators, which depend 
on discrete variable k, for the operators, depending on 
continuous variable z, that is, 



(21) 



Thus, we obtain, taking also into account the relation- 
ships 

2 d 2 &^(z) 



a 



z,± 



the equation, which, like to (|18p . in compact vector form 



2,± 



(22) 



at 

where 



,r(z) x 7e £ 



2a 2 J„ 



<?(z) x V 2 ct(z) , (23) 



E = E 1 e wt e-+E 1 e~ 
= E+ ( 



-LOq 



(24) 



E z e 2 



The structure of vector E clarifies its physical meaning. 
Two components E + , E~ are right- and left-rotatory 
components of oscillating external electric field, third 
component E z is intracrystalline electric field, which pro- 
duces two level energy splitting for each of the elementary 
unit of a chain system with value, equal to Hluq. When 
using the analogy with magnetic resonance the following 
relationship can be entered 



(25) 



where f3 E is electrical analogue of Bohr magneton, g E 
is electrical analogue of magnetic g -tensor, which is as- 
sumed for the simplicity to be isotropic. In other words, 
by means of relationship [25] the correspondence between 
an unknown intracrystalline electric field Eq and ob- 
served frequency ujq is set up. Further, we take into 
consideration, that physical sense of operators <y(z) in 
continuous limit is remained, that is, for each point of z 
the components a ± (z), o z [z) are satisfying to algebra, 
which is isomorphic to algebra of the set of spin compo- 
nents. Then by means of relationship, which has mathe- 
matically the same form for both the types of dynamical 
systems studied S(z) ~ hcr ( z ) ; ^he equation of the motion 
for operators of magnetic and electric spin moments are 
obtained. So, for example, the equation of the motion 
for electric spin moment operator is: 



dS(z 
dt 



S(z) x 1e E 



4« 2 J E 

h 2 



S{z) x V 2 S(z) (26) 



III. 



L-L EQUATION SOLUTION AND 
DISCUSSION 



We see, that transition operator method, being apply- 
ing to optical transitions, leads unambiguously to con- 
clusion on existence of electric spin moment, that is to 



the same conclusion, which was done by Dirac in 1928 
in its famous paper although Dirac himself was in 
some doubt about physical and even mathematical cor- 
rectness of given conclusion. Transition operator method 
relieves all doubts in mathematical correctness of Dirac' 
conclusion (although method does not give any indica- 
tions about imaginary or real character for given quan- 
tity). Given aspect was discussed in recent paper [lol ]. 
where experimental confirmation for electric spin mo- 
ment discovery was represented. 

Equation (|2"6")) gives for the case J = quantum- 
mechanical optical analogue of classical Landau-Lifshitz 
equation in continuous limit. In fact equation (|26[) and 
its magnetic analogue (if replace electrical characteris- 
tics by corresponding magnetic characteristics by means 
of eq. (fl~3"|0 are operator equations, which argue the cor- 
rectness of eqs.(Q] [2]), that is, the physical correctness of 
gyroscopic model for description of spectroscopic transi- 
tions and transitional phenomena for both optical and 
magnetic resonance spectroscopy. If J ^ equation ([26]) 
gives quantum-mechanical operator analogue of the equa- 
tion for observables, which was introduced by Kittel for 
ferromagnetic spin wave resonance (FM SWR) descrip- 
tion 0]. The equation (f26j) (if put aside the operator 
symbols and to replace electrical characteristics by cor- 
responding magnetic characteristics) and the equation, 
introduced by Kittel, are coinciding mathematically to 
factor two in the second term. This difference is like to 
well known difference of gyromagnetic ratios for orbital 
and spin angular moments. Given comparison means, 
that equation (f2"6"| is the equation for optical analogue of 
ferromagnetic spin wave resonance. In other words, we 
have obtained the equation, which predicts new physi- 
cal phenomenon - ferroelectric spin wave resonance (FE 
SWR). experimentally and reported in [111 ]. 

We can also conclude, that the spin value can be reg- 
istered in solids directly by optical methods. Really, it 
will be possible in the cases, when the responce of dy- 
namical system studied can be obtained by registration 
schema, which takes into account the imaginary nature of 
electric spin moment. Given schema has to ensure -| hy- 
perbolic phase rotation in the functional space, which is 
isomorphic to Minkowski space like to microwave ? phase 
rotation by detection of the signals, being proportional 
to real and imaginary parts of complex magnetic suscep- 
tibiliy in electron spin resonance experiments. On the 
other hand spin registration can be provided automati- 
cally, if in Hamiltonian there is item with product of the 
components of electric spin moment like to Hamiltonian 
(fT2"]) . Consequently, transition operator method is the- 
oretical substantiation for experimental methods of spin 
detection. Thus along with magnetic resonance methods 
we can detect a spin value of particles, quasiparticles, 
impurities or other centers in solids by optical methods, 
for example by study of transitional optical analogues of 
magnetic resonance and especially simply by new pre- 
dicted phenomenon FE SWR. 

It should also be noted, that above considered theo- 
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retical description of FE SWR allow predict the differ- 
ence in splitting constants which characterize FE SWR 
by its experimental detection with using of one-photon 
methods like to IR-absorption or IR-reflection and with 
using of two-photon methods like to Raman scattering. 
It is evident, that equations ([23]) . (|26|) can immediately 
be used for single transition methods, for instance, for 
IR-absorption. By Raman scattering we have two sub- 
sequent transitions. Then operator <j(z), which char- 
acterizes the transition dynamics by Raman scattering 
process, has to be consisting of two components <?i(z) 
and a 2 {z) characterizing both the transitions, taken sep- 
arately, that is &{z) = o*x(z) + 0*2(2). Consequently, the 
equation for transition dynamics of second component, 
which will determine experimentally observed FE SWR- 
spectrum, is 



da 2 (z) 
dt 



a 2 (z) x "f E E 
2a 2 J 



2a 2 J 



a 2 (z) x V 2 a 2 (z) 



0*2(2) X V 2 CTi(z) 



(27) 



The second and the third items in (|III[) are practically 
equal to each other, since we are dealing with interact- 
ing electric dipole moments of the same chain, that is, 
V 2 ai(z) and V 2 o*2(z) have almost equal values. Then 
we obtain, that the value of splitting constant by Raman 
scattering detection of FE SWR in the same sample is al- 
most double in comparison with that one by IR detection 
of FE SWR. The observation of doubling in the splitting 
constant by Raman FE SWR-studies is additional direct 
argument in FE SWR identification. 

FE SWR was discovered experimentally and reported 
in [111 ]. To compare with experiment the equation (112[) 
and consequently final equation (|26[) have to be modified 
by taking into consideration the relaxation processes and 
from operator variables we have to proceed to observ- 
ables. Modification was done in [ll[ fenomenologically 
by adding the term, which is similar to corresponding 
term in Bloch equations. If J 7^ equation ([21?]) gives 
quantum-mechanical operator analogue of the equation 
for observables, which was introduced by Kittel for ferro- 
magnetic spin wave resonance (FM SWR) description Q ■ 
It can be seen, that the coincidence of equation (|26[) and 
Kittel equation Q will take place, if put aside the oper- 
ator symbols and to replace electrical characteristics by 
corresponding magnetic characteristics. Consequently, to 
proceed to observables it is sufficient to put aside the op- 
erator symbols. It was suggested in Then, for de- 
scription of optical transition dynamics the equation will 
be 



4a 2 J E 



dS(z) 
dt 

S(z) x W 2 S(z) 



S{z) x lE E 
S{z) - S (z) 



(28) 



where Sq(z) is equilibrium value of electrical spin mo- 
ment vector function, r is relaxation time. It was also 



suggested in [Tl| that experinental data on FE SWR, 
which was detected by three various optical methods - 
by IR absorption, by IR reflection and by Raman scat- 
tering - will be well described by linearized equation, cor- 
responding to equation (]26p . The linearized equation in 
11] is considered. It was obtained under the assump- 
tion, that the values of oscillating external electric field 
components E x , E y in (f2"8"|) are in experiment greatly less 
in comparison with the value of intracrystalline electric 
field component E z , under analogous assumption rela- 
tively the components of total electrical spin moment and 
under additional assumption, that equilibrium distribu- 
tion of So(z) along the chain is homogeneous. All the as- 
sumptions are entirely correct for IR measurements. The 
linearized equation was solved and the relationships for a 
shape and amplitudes of resonance modes and dispersion 
law were obtained. They are 



»7e St 2 E x 



[(<* 



■j)- 



mi 



1, 3, 5,.. 



0. 



n = 2, 4, 6, . 



vo — 2ln , 



(29) 



(30) 



where n € N including zero, v n is a frequency of n-th 
mode, 21 is a material parameter (21 = 27r ° lL l^ n 2 > 0). 
Here Re a n is proportional to absorption signal, Im a n is 
proportional to dispersion signal. It was found the fol- 
lowig. 1) Dispersion law ([30]) is held true both by IR- 
and RS-detection of FE SWR. 2) The excitation of the 
only uneven modes by IR- detection (by which the exper- 
imental conditions were corresponding to applicability of 
linearized equation) takes place in accordance with (|29|) . 
3) Inversely proportional dependence at resonance of the 
amplitudes of modes on mode number n in accordance 
with (|29p is also held true, however also the only by IR- 
detection of FE SWR. 4) Splitting of Raman active vi- 
bration modes is characterized by value of parameter 21, 
being approximately by factor 2 greater, than parame- 
ter 21, which characterizes IR FESWR spectra (by the 
frequencies of zero modes reduced by means of linear ap- 
proximation procedure to the same value). The fourth 
result agrees well with the above discussed results, pre- 
liminary reported in Q. The agreement of theoretical 
and experimental results, which is rather well, allow to 
obtain new substantial information on carbynoid struc- 
ture. Since Hamiltonian (TT2"]) is true the only for the 
centers with the spin of i it means that spin-Peierls soli- 
tons [13], [HI both optical and magneticai, which are 
formed in carbynoid chains have spin S = \. The obser- 
vation of two types of FM SWR extitations [l3| and two 
types of FE SWR extitations with slightly different split- 
ting parameters 21 ff and 21 £ indicate unambiguously that 
spin-Peierls solitons are produced separately in tt x — and 
^-subsystems of carbynoid electronic 7r-system. More- 
over we can also conclude, that tt x - and ^-subsystems 
of carbynoids are inequivalent. It seems to be essential, 
that given inequivalency was observed after 12 Y stor- 
age. The only one set of FM SWR has been found in the 



samples, which were studied immediately after their pro- 
duction ([l2j , [HI), indicating on full equivalency of tt x - 
and 7r y -subsystems of starting carbynoids. Naturally, all 
electronic properties, wich determined mainly by tt x - and 
•^-subsystems are also changed with storage time. How- 
ever, clear above presented elucidation of the cause for 
the changes of electronic properties gives key to produc- 
tion of storage stable carbynoids. The values of splitting 
parameters 21 B and 2l H allow to find the ratio Je/Jh 
of exchange constants. The values of splitting parameter 
% H of FM SWR which was observed strictly in the same 
sample, that is on the same chains, were reported earlier 
in pj. The range of the ratio J e /Jh is (1.2 - 1.6)10 4 . 
Given result seems to be direct proof, that the function, 
which is invariant under gauge transformations is two 
component, that is complex- valued function. In other 
words, the complex charge corresponds to presence of 
exchange interacting solitons with two independent ex- 
change constants. We can evaluate the ratio of imagine 
en = g t o real eg = e components of complex charge to 
be § ~ ^JeIJ^ « (1.1 - 1.3)10 2 . 



da+(z,t) 
dt 



iuj a + (z, t) + ifl R e luJt a z (z, t) 



(31) 



+ 



2ia 2 J\^, „d 2 a z (z,t) d 2 a+(z,t) 



a+(z,t)- 



dz 2 



a z (z,t)- 



dz 2 



da-(z,t) 
dt 

2ia 2 J 



= -iu} Q a~ (z,t)- in R e- luJt a z (z, t) 
- .. 9VM ,,d 2 a-(z,t) 



(32) 



dz 2 



dz 2 



da z {z,t) 

dt 
Aia 2 J 



= 2 



itt R (cr+ (z, t)e~ wt -a~(z, t)e iut ) ^ 



. ,d 2 a+(z,t) ,, d 2 a~(z,t) 
a-(z,t) -\^-a+{z,t) 



dz 2 



Let proceed to new variables 



dz 2 



dS'*(z,t) 
dt 

Aia 2 J 



= 2iSl R [S + {z,t)-S-(z,t)] 
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(37) 



dz 2 



dz 2 



where A = loq — to. 

Passing on to new basis e x , e v , e z and designating 



S x (z,t) 



2a 2 J 
.2a 2 J , 



(S+(z,t) + S-(z,t)), 



S y (z,t) = i — (S-(z,t) - S+(z,t)), 



S z (z,t) 
we obtain 

as* CM) = 

dt 



4a 2 J 



r, 2 



S' z (z,t), 



ASy{ Z ,t) 



dz 2 



dz 2 



dSy(z,t) 

dt 



AS x (z,t) + n R S z {z,t) 



^ d 2 S z (z,t) ,,d 2 S x (z,t) 
S x \z,t) — S z {z,t) 



dS z (z,t) 
dt 



dz 2 



= -Q R Sy(Z,t) 



dz 2 



^ d 2 S x (z,t) ^ d 2 S y (z,t) o 



dz 2 



dz 2 



or in vector form 
dS(z,t) 



dt 



S(z,t) x 



d 2 S(z,t) 
dz 2 



S(z,t) x n , (42) 



(38) 



(39) 



(40) 



(41) 



S + (z,t) = -a+(z,t)e~^ 4 
S-(z,t) = ^a-(z,t)e^, 
S' z (z,t) = ^* z (z,t), 



then we will have 
dS+(z,t) 



dt 

Ma 2 J 



dS~(z,t) 
dt 

4ia 2 J 



iAS + (z,t)+m R S' z (z 7 t) 

d 2 s ,z {z,t) 



(34) 



(35) 



S + (z,t) 



dz 



S' z {z,t) 



d 2 S+(z,t) 



dz 



r, 2 



= -iAS~(z,t) - ifl R S' z (z,t) 



(36) 



dz 2 



dz 2 



S(z,t) = S x (z,t)e x + S y {z,t)e y + S z {z,t)e z , (43) 



SI = Q R e x — Ae* z . 



(44) 



The equation (??) has exact analytical n-soliton solution 
[LH ]. In particular, at n = 1, the solution is 

S x (z, t) = sin 9 cos <p, S v {z, t) = sin 9 sin ip, S z (z, t) = cos#, 

(45) 

where 



cosfl = 1 - 20 2 chT 2 \j3yfu{z -vt- z )], 
W = 4i?e[£],^ = 4|£| 2 , 

ip ~ ipo + tot + v{z — vt) /2 
+ arctg{f3(l - (3 2 )- 1/2 th[f3y/^(z - vt - z )}}, 



(46) 



(47) 
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= z = Zri{M/(2Im[£])]}/(2/m[£]), 
Im [£,]/\£.\, ( Po = arg[m}. 



(48) 



Here £ is eigenvalue, m is asymptotic characteris- 
tic of eigenfunction of differential operator L(z,t) = 

(S(z,t)ap)d z , op is Pauli matrix vector. The soliton 
is characterised by four real parameters: velocity v, fre- 
quency u>, phase ipo and coordinate of the center Zq. 

It should be noted, that cosideration of both variants of 
solutions, that is exact solution and solution of linearized 
equation, seems to be necessary. It follows from general 
virial theorem, which recently was proven by Buniy and 
Kephart. For any classic field, independent of a partial 
model, the theorem establishes, in particular, the lower 
energy bounds for massive classical solitons or lumps [l5| . 
It means, that soliton solution will be applicable, when 
EM-field energy exceeds some threshold. 

Let us remember that exact solution and solution of 
linearized equation were obtained by means of proceed- 
ing to observables in operator equation in s simple way 
- by change of operators into corresponding observables. 
Agreement with Kittel theory and, the essentials, agree- 
ment with experimental data justify given suggestion, 
however mathematical proof is desirable. On second 
hand EM-field is relativistic object. It means, that the 
equations, describing the processes of interaction with 
EM-field, strictly speaking, have to be Lorentz-invariant. 
To solve given tasks we use further quaternion formalism. 
When applying to quantum mechanics quaternion for- 
malism realizes its generalization, consisting in replacing 
of standard Hilbert space over the field of complex num- 
bers into Hilbert space over the ring of quaternions. More 
strictly, biquaternions are used. We use the quaternion 
basis {ei},i = 0,3 with algebraic operations between el- 
ements, satisfying to relationships 

^ijk&k ~t~ &ij &o ) &i ^ij ' J> ^ ^' ^> 

. . (49) 

where Eijk is completely antisymmetric Levi-Chivita 3- 
tensor. 

Then we use the equation completed like to pi))) 
by item, taking phenomenologically the relaxation pro- 
cesses 



da(z,t) 
8t 

2a 2 J ' 



a(z,t) x F(z,t)E 
<r(z,t) x V 2 a(z,t) 



a(z,t) - a (z) 



(50) 



where F — j E E, E is unit operator. Quaternion algebra 
and algebra of Pauli matrix, completed in the manner 
above described, are isomorphic, see for example [ijj]. It 
means, that algebra of transions operators is also isomor- 
phic to algebra of complete Pauli matrix set. It is rea- 
sonable to introduce the quaternion operators by means 
of mappings E — > e ,a(z,t) — > ie, e — (ei,e2,e3) with 
conditions V/(cc) the relationships e^f{x) = ekfk{x) is 
taking place. Here <I>(x) = f (x)e + fi(x)e{ is arbitrary 
quaternion function of quaternion argument 



x = ict, that is, components of Xi G 1 R4. Let us define 
now wave function $(x) of the chain system interactig 
with EM-field and described by Hamiltonian (fT3i|) to be 
vector of the state in Hilbert space over quaternion ring. 
In other words, it is quaternion function of quaternion 
argument and can be represented in a form 
$>(x) = $ (x)e +Phi l (x)e t . 

It reazonable to suggest, that $ (x) = 0. Then taking 
into account replacing t — > ict = 
be rewritten in the form 



equation (|50|) can 



dx 
2a 2 J, 



£\®(x)\ x F(x)e 



x V 2 (^$(a;)|) 



(51) 

g|*(aO|-g|* (aO| 



where |$ (a;)| is |$(x)| at t = t , F(x) is vector quater- 
nion delivered in conformity with vector F{z, t). It is is 
evident, that e|$(a;)| = $(x). Then we have 



d$(x) 
dx 
2a 2 J, 



1 r 

ic 



lie 



<f>(x) x F(x) 
$(x) x V 2 $(x) 



$(x)-$ (x) 



(52) 



The equation obtained is Lorenz invariant equation for 
vector of the state <&(x) of the system EM-field + chain. 
Quaternion vector of the state <f>(x) is transformed ac- 
cording to vector representation of Lorenz group, which 
is $'(x) = L<fr(x)L*. Spin vector S(x) is transformed 
under the same vector representation of Lorenz group. 
So it means, on the one hand, that ([52]) will held true, 
if <i>(x) replace into S(x), and we obtain therefore the 
equation for observables. On the other hand, S(x) can 
be considered to be quaternion vector of the state of the 
system studied. 



IV. CONCLUSIONS 

Therefore, quantum-mechanical analogue of Landau- 
Lifshitz equation has been derived with clear physi- 
cal meaning of the quantities for both radio- and opti- 
cal spectroscopy. It has been established that Landau- 
Lifshitz equation is fundamental physical equation un- 
derlying the dynamics of spectroscopic transitions and 
transitional phenomena. The solution of both linearized 
and starting nonlinear equation is presented. New phe- 
nomenon - ferroelectric spin wave resonance and its some 
spectroscopic peculiarities are described. The agreement 
of theoretical and experimental results, which is rather 
well, allows to obtain new substantial information on 
chain structure and the properties of chain elementary 
units. For example, for the first time from pure optical 
measurements the value of spin 5 = 1/2 for optically 
active centers - spin-Peierls solitons has been obtained. 
The chains of spin-Peirls solitons are exchange coupled, 
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the interaction is characterized by complex- valued charge 
function. The ratio of imagine to real components of 
complex charge is evaluated to be | w (1.1 — 1.3)10 2 . 

Lorenz invariant equation for vector of the state of the 
system EM-ficld + chain has been obtained and solved. 
It is shown, that spin vector can be considered to be 
quaternion vector of the state of the system studied. 
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